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In the theory of open quantum systems interaction is a fundamental concepts in the review of the dynamics 
of open quantum systems. Correlation, both classical and quantum one, is generated due to interaction between 
system and environment. Here, we recall the quantity which well known as total entropy production. Appear¬ 
ance of total entropy production is due to the entanglement production between system an environment. In this 
work, we discuss about the role of the total entropy production for detecting non-Markovianity. By utilizing the 
relation between total entropy production and total correlation between subsystems, one can see a temporary 
decrease of total entropy production is a signature of non-Markovianity. 
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I. INTRODUCTION 

The study of open quantum systems from different perspec¬ 
tives has been the interest of many researchers [1-3]. Infor¬ 
mation flow between the system and its surroundings is one of 
the main characteristics of the open quantum systems. Indeed, 
the flow of information is due to the interaction of the system 
with its surrounding environment which its direction depends 
on several factors such as the amount of coupling between 
system and environment. From the standpoint of the mem¬ 
ory effects, the quantum dynamical processes is divided into 
two categories, namely, Markovian (memoryless) and non- 
Markovian (with memory) dynamical maps. For Markovian 
dynamical maps (memoryless process) the coupling between 
system and environment is weak and information flow from 
system to environment continuously. If the coupling between 
system and environment is strong the memory effects will be 
apparent and the future states of the system depend on its past 
which it is the result of the back flow of information from 
environment to the system. There is much to learn from the 
presence of the past in the future, because of this the study 
of open quantum system in non-Markovian regime, detecting 
and quantifying that by several distinct criteria , are the major 
goals of the new researches [5-9]. One can see the completely 
positive dynamical map A is markovian, when it forms a one 
parameter semi-group with generator C in Lindblad form [1] 

Co=-i[H, o] + 5>n(2V„oV2 _ {V n V n ,o}), (1) 

n 

where 7 n > 0 for every n, H is Hermitian operator and 
V n ’s called Lindblad operators acting on the system’s Hilbert 
space. In the case of time dependent, C t is referred to time de¬ 
pendent Markovian evolution when 7 n (t) > 0 for every t and 
n, note that here C t does not lead to one parameter semi-group 
of dynamical maps. Due to this feature violation of semi¬ 
group property is not sufficient for dynamical map to be non- 
Markovian, extra property well known as divisibility, must be 


violated [4]. Divisibility directly corresponds to Markovian 
dynamics and non-Markovianity can be quantified as the de¬ 
gree of deviation from divisibility of dynamical maps, based 
on this feature, Rivas et al. proposed a measure well known as 
RHP [5]. Another phenomena which is used to quantify non- 
Markovianity is the reduction of distinguishability of quan¬ 
tum states under completely positive trace preserving maps, 
which is associated to a loss of information about the quan¬ 
tum system. Utilizing this fact, Breuer, Laine and Piilo de¬ 
fined a measure which is called BLP [7]. A series of measures 
have been introduced in the light of the behaviour of correla¬ 
tions, both classical and quantum one, under Markovian and 
non-Markovian dynamical maps [5, 8, 10, 11]. For example, 
one can refers to the criteria has been proposed by Rivas et 
al. which is based on the entanglement between a system and 
an isolated ancilla [5]. Fanchini et al. By using accessible 
information give an interesting interpretation for this criteria 
[11]. In the same manner Luo, Fu, and Song, use mutual infor¬ 
mation to detect non-Markovianity [8], their measure have an 
information interpretation in the context of information loss, 
which is introduced by Haseli et al. [12]. As we know creat¬ 
ing correlation, both classical or quantum one, is one result of 
interaction between system and environment, which is caused 
disorders in total system. During the interaction total entropy 
changes and amount of this change is always positive and it is 
called total entropy production ( TEP ) as a measure of disor¬ 
ders [13, 14]. Appearance of TEP is due to the entanglement 
production between system an environment. In this work ac¬ 
cording to decoherence model, where an isolated system S is 
coupled to a measurement apparatus A, which in turn inter¬ 
acts with an environment £, we want to utilize TEP quan¬ 
tity in orther to introduce a measure for non-Markovianity. 
First, we obtain an lower bound for total correlation changes 
A I(S: A ), interestingly this lower bound in associated with 
TEP by an inequality. One can see, if total entropy produc¬ 
tion decrease monotonically increasing in total correlation is 
inevitable which is a signature of non-Markovianity. This pa¬ 
per is organized as folio. In Sec. II we introduce the defini¬ 
tions of the entropy exchange. In section III, by making use of 
the connection between total entropy production and quantum 
mutual information, we introduce the definition of the related 
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non-Markovianity measure. In section IV, some examples in 
dynamical model are provided in order to examine the mea¬ 
sure. Finally in section V, we summarize our discussion and 
results. 


II. ENTROPY EXCHANGE 



Entropy exchange was introduced by Schumacher [13] and 
Lloyd [16]. The entropy exchange of quantum operation A, 
with input p is defined to be 


FIG. 1: We consider an initially pure environment 8, and an en¬ 
tangled pure state SA. As the system S evolves free of any direct 
interaction, the apparatus is interacting with the environment 8. 


S e (p,A) = S(A[p]) = S(/), (2) 

where p s is the state of the environment after the operation 
A. Note that the state of the environment initially is assumed 
pure. If A [p\ = JA FipFi then an appropriate form for en¬ 
tropy exchange by introducing a matrix W with following el¬ 
ements 


Wij = 


TrjFipF)) 

Tr(A[p\) ’ 


( 3 ) 


can be written as 


here Jy{p xy ) = S{p x ) - min {n y } Y Ji p i S(p x \ y ) is 

the classical correlation between X and y if an ob¬ 
server doing measurement on y and Dy(p xy ) = 
— S{p x \ y ) denotes quantum dis¬ 
cord with respect to measurement on y, where S(p) = 
— Tr(p\og 2 p) and S(p x \ y ) = S(p xy ) — S(p y ) is the von 
Neumann and conditional von Neumann entropy respectively. 
Using the definition of the quantum discord in our decoher¬ 
ence program the amount of change in the quantum correla¬ 
tion can be derived as 

A D(S: A) = D A (p SA ) - D A (p SA ) (7) 


S e (p, A) = S(W) = - Tr(W log 2 W). (4) 


= [S(p s I' 4 ) - S(p s !■*)]+ 


III. ENTROPY PRODUCTION 

In this section, we introduce a new method of quantifying 
non-Markovianity through the total change in the entropy of 
the total system SA and environment^. Our method is estab¬ 
lished on the decoherence program, where a quantum system 
S is cou- pled to a measurement apparatus A, which in turn 
directly in- teracts with an environment 8. One can consider 
a quantum system S that is initially correlated with the appa¬ 
ratus A and the state SA is in product state with environment 
8. Initially, the state SA , total state SA8 and the state of the 
environment 8 is assumed pure. The environment only affects 
the state of the apparatus A. As a result of the interaction, 
there emerges an amount of correlation among the individual 
parts of the closed tripartite system SA8 , and thus the envi¬ 
ronment 8 acquires information about the system S by means 
of the interaction with the apparatus A. This setting is graph¬ 
ically sketched in Fig. 1, where the system S evolves trivially 
while the apparatus A is in a direct unitary interaction with 
the environment E . The final state of the composite bipartite 
system SA is given by 

p SA = Tr £ [(I®U A£ )p SA£ (I®U A£ Y} = [I®A}p AS , (5) 

A is the general completely positive trace preserving map. Let 
us now describe our strategy to derive our witness to detect 
non-Markovianity. For a bipartite system p xy quantum mu¬ 
tual information is defined as 

( 6 ) 


+ min {n A } Y,PiS(p SlA ) ~ min {n A } ypjS(p s ^ A ). 

i i 

We introduce quantity that will be very useful, namely, to¬ 
tal entropy production(TEP) [13-15]. The difference between 
final total entropy Sp = S(p SA ) + S(p^) and initial total 
entropy Sp = S(p SA ) + S(p £ ) is defined as total entropy 
production TEP 

A S P = Sp-S P = A Ssa + AS*. ( 8 ) 

Note that (TEP) is always positive. By using the initial and 
final total entropy and substituting Eq .8 in to Eq.7 we have 

A D(S: A) = -AS P - [J{p slA ) - J(p SlA )}~ (9) 

- [S(p A£ ) - S(p A ) - S(p £ )}. 

According to the definition of mutual information in Eq. 6 , one 
can rewrite Eq.9 as 

A I(S: A) = -A S P - [S(p A£ ) - S(p A ) - S(p £ )], (10) 

where A I(S: A) = I - I = A D(S: A) + A J(S: A) is the 
total change in mutual information. Using the fact that A 
and 8 initially are in product state and subadditivity of the 
von Neumann entropy and since p AE evolve unitarily we have 
S(p A£ ) < S(p A ) + S(p £ ) andEq.lOleadto 

A I(S:A) > -A S P , 


I(p Xy ) = Dy(p Xy ) + Jy(p Xy ), 


(ID 
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for very small timescales Eq. 11 can be rewritten as 


IV. SOME DYNAMICAL EXAMPLES 


dl dS P 
dt ~ dt 


( 12 ) 


A. Pure dephasing model 


Recalling that the LFS measure of non-Markovianity is based 
on the rate of change of the quantum mutual information 
shared by the system S and the apparatus A. In particular, the 
LFS measure captures the non-Markovian behaviour through 
a temporary increase of the mutual information of the bipartite 
system SA, from Eq.12 this satisfy if 


dS P 

dt 


< 0 . 


(13) 


It captures the non-Markovian behavior through a temporary 
decrease of the (TEP) of the tripartite SAS. Mathematically, 
the measure can be written as 


A/sp = 


max 



dSp 

dt 


dt , 


(14) 


where the maximization is evaluated over all possible pure 
initial states of the bipartite system SA. Due to the fact that 
the initial state of tripartite SAS is pure and remain pure be¬ 
cause of unitary evolution, the von Neumann entropy of each 
subsystem SA or S represents the amount of entanglement 
between them, hence Sp = 2 S(p s ), which is equal to twice 
amount of entanglement between SA and £ at each instant 
of the process. Physically, one can say that, if the amount of 
entanglement between SA and £ decreases then we confront 
with the non-Markovian processes. As a result of the Eq. 12, 
the process is non-Markovian iff 


dSjpf) _ dSjp |S - 4 ) 
dt dt 


(15) 


One can employ spin-boson model to construct a general 
pure dephasing model [1-3, 17]. We consider a two level 
quantum system as an apparatus A, which linearly interact¬ 
ing with an environment £, so that the total hamiltonian is 

H = + ^2,uJib\bi + cr z {gib\ + g*bi). (16) 

i i 

we will focus on the case of initial pure state of the bath. In 
order to continue the procedure , we have to know the spectral 
density J(cc) of the bath. Here we consider the ohmic spectral 
density as follow 


J(^) = A=■ ( 17 > 

UJ C 

uj c is cut-off frequency of the spectrum and s is bath param¬ 
eter. According to this parameter, the bath change from sub- 
ohmic s < 1 to Ohmic s = 1 and super-Ohmic s > 1 [18]. 
Under these conditions, the reduced dynamic is calculated to 
give 


o(t) _ ( Pn Pi2l(t)\ 

P[t) ~ Ul7C0 P22 ) 

where the dephasing parameter 7 (t) is 

7 (t) = exp[- / ??(r)dr], 


where dephasing rate rj(r) is 


(18) 


(19) 


where S(p s ) is the entropy exchange. 



FIG. 2: (Color online) Total entropy production as a function of 
scaled time 70 1 (red line). The inset displays the behaviour of total 
entropy production rate dASp/dt (blue line) under pure dephasing 
channel with super Ohmic reservoir spectrum by s = 4 


rj(r) = u) c [l + (cc c t) 2 ] 2T(s) sin[sarctan(o; c r)], (20) 

with T(s) the Euler gamma function. In Fig. 2, we can see the 
plot of the total entropy production rate dAS p /dt ( TEPR ) 
(red line), the inset displays the behaviour of TEP (blue line) 
as a function of scaled time Xt . in this model. As one can seen, 
the negative value for toral entropy production rate TEPR is 
appeared for this model by considering special condition on 
Ohmicity parameter s. Note that in this model dynamic does 
not admit non-Markovianity for s < 1. 

B. Amplitude damping 

Here we consider the apparatus A as a two-level quantum 
system which interacts with zero temperature environment. 
The dynamics is governed by the following interaction Hamil¬ 
tonian 

H = u 0 cj + (j-+'Y^io k a\a k +'Y^(g k (T + a k +gl(T-a\) (21) 

k k 

where cco is the transition frequency of the apparatus A and 
a± denote the raising and lowering operators of the apparatus 
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A. dk and a\ are the annihilation and creation operators of 
the environment £ , respectively, with the frequencies ujk • The 
spectral density of the environment has the Lorentzian form 

J{uo) = 7 0 A 2 /27t[(u;o - ^) 2 + A 2 ], (22) 

where the A is the spectral width of the coupling. A is related 
to the correlation time of the environment by ~ 1 / A 
and 70 is related to the time scale tr by tr « I/ 70 , where 
tr is the time scale that in which the system change. Accord¬ 
ing to these considerations We can define the dynamics of the 
apparatus via the following master equation as 

^p(t) = 7 (t) (cr_p(t)cr + - t{ a+< j_,p(t)}) , (23) 



where 7 (t) is the time-dependent decay rate and is given by 
27 0 Asinh (dt/2) 


7 (t) = 


d cosh (dt/2) + A sinh (dt/ 2 ) ’ 


(24) 


where d = ^/A 2 — 270 A. Thus we can define the dynamics of 
the A via the Kraus representation as 


P(t) = A t , p=Y, Ki(t)pKi(tf, (25) 

km 1 

the corresponding Kraus operators is given by 
Kiit) = 


1 0 
0 G(t) 


) > k 2 (t) = ~WW) ; 


where the function G(t) has the following form 


G(t) = e 


= P~ xt / 2 


, dt > A , dt , 

cosh( —) + ^ smh (y) 


(26) 


(27) 


FIG. 3: (Color online ) Total entropy production as a function of 
time (red line) and the inset shows the total entropy production rate 
dAS p /dt (blue line) for the exact amplitude damping model with 
Lorentzian spectral density for reservoir with A /70 = 0.05 



FIG. 4: (Color online ) Total entropy production as a function of 
time (red line) and the inset shows the total entropy production rate 
dASp/dt dynamic (blue line) for the generalized amplitude damping 
model with uj = 5. 


As can be seen from Fig. 3, by choosing the value A /70 = 
0.05, total entropy production rate TEPR in some intervals of 
time takes the negative value and we have the temporary loss 
of total entropy production TEP, it shows that in this situation 
the dynamics is non-Markovian which is consistent with the 
results from other measures [5, 7-9, 11, 12]. 


C. Generalized amplitude damping 

In this case, we consider a two level system as an apparatus 
A which interacts with an environment at finite temperature. 
A generalized amplitude damping channel describes the re¬ 
laxation due to coupling of system to their surrounding with 
temperature often much higher than the system [19]. For sin¬ 


gle qubit systems it is defined by following kraus operators 

EiM = dm (/ y^y), 

E 2 {i) = y/m (jj ^ “ ^) , 

E 3 (t) = y/l^W) (Tjp) 0) , 

E 4 (t) = yrrpw jj), (28) 

whare for all time t we have Ej(t)Ei(t) and 

P(£), q(t) G [0,1]. For simplicity we choose the parameters 
as P(t) = cos 2 cct and q(t) = e _t , where uo is a real num- 
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ber. In Fig. 4, one can see for some different time intervals the 
TEPR is negative, hence TEP decrease and the dynamics is 
non-Markovian [12]. 

V. RESULTS AND DISCUSSION 

In this paper, our discussion relied on the decoherence 
model, where the system S was coupled to the apparatus A 
and initially the composite system SA was correlated and in 
pure state. The environment 8 directly just had an interaction 
with apparatus A and the system S remain intact. We began 
our strategy with total correlation between system S and ap¬ 
paratus A and continued by the difference between final and 


initial quantum discord between them. Finally, after straight¬ 
forward calculation we obtained the inequality between the 
time derivative of mutual information and entropy production 
which by making use of the Luo measure [8] we defined the 
new criteria in order to detect non-Markovianity based on en¬ 
tropy production. Here, we can give a physical interpretation 
to this criteria. Because of the interaction between environ¬ 
ment 8 and apparatus A the decoherence will be occur and 
we have the entanglement production between <5^4 and 8. If 
the amount of the produced entanglement between SA and 
8 decreases then we confront with the non-Markovian pro¬ 
cesses. 
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